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Abstract 

We examine the consequences of leptoquark explanation of HERA anomaly in 

the context of R parity conserving supersymmetric gauge unified theory with 

the gauge unification scale at ~ 10^^ GeV. We pointed out the difficulty of 

constructing a grandunified theory. However gauge unification is still possible 

at ~ 10^^ GeV when additional multiplets are introduced. We determine 

the mass spectrum of these additional fields (fermions and scalars) in gauge 

mediated and supergravity scenarios. Unique signatures and mass bounds are 

discussed. 
PACS numbers: 
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Recent observation of high anamoly in e~^p scattering seen by HI |T| and ZEUS 
collaboration at HERA, if confirmed by further experimentation, would be a clear indica- 
tion of physics beyond the Standard Model (SM). Explanation of this anomaly have been 
advanced in terms of contact interactions generated by physics beyond the SM , in terms 
of s-channel leptoquark production |^ or in terms of modification of the parton distribution 
function DH- Contact interactions have to satisfy various constraints e.g. atomic parity vio- 
lation, LEP data, CCFR data etc. Although it is possible to construct such models, these 
typically predict similar anomalous events in e~p data and excess e~^e~ events in Drell-Yan 
process at high in FERMILAB, evidence for either is lacking at present. Leptoquark 
explanation has been pursued in the context of R-parity conserving or R-parity violating 
supersymmetric theory. These states occur as resonances at a mass of ~ 200-220 GeV. In 
this paper we shall pursue the leptoquark explanation of HERA data, but in the context of 
R-parity conserving supersymmetric gauge unification. We prefer to maintain R invariance 
because this allows the lightest SUSY particle to be a candidate for dark matter. Further, in 
R violating scenarios there are numerous arbitrary coupling e.g 27 of X'{\[ji^LiQjD'^), 9 
of \{\ijkLiLjE'^) and 9 of \"{X'-jf.UfDjDj(), where Lj, Qi are the lepton and quark doublets 
and E'^, W^, are the lepton, uptype and down type singlets respectively; and also A' and 
A" can not appear together in the same superpotential, since that gives rise to rapid proton 
decay. In this present paper we shall demand that any additional leptoquark states pos- 
tulated in the SU{?>) x SU{2) x U{1) context be part of supersymmetric gauge unification 
picture so as to be theoretically attractive. Since the leptoquark mass is 200-220 GeV, we 
expect that there are other fields around this mass in order to maintain gauge unification. 

We first consider gauge unification in the context of a grandunified theory like SU(5) 
0. Lowest dimensional representation of SU(5) that we can have leptoquarks in is the 
10 dimensional representation. However, because 10 dimensional representation is formed 
from the antisymmetric combination of 5 x 5, its SU(5) invariant coupling 5i5jl0k is an- 
tisymmetric in the generation indices i and j. Thus the leptoquark can not couple to the 
first generation quarks and lepton as needed to explain the HERA events. One then has 
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to consider a higher dimensional representation hke 15, which is symmetric in the product 
5x5. To remove gauge anomahes one also requires 15 representation, which also makes 
the representation vector like. However 15 representation has diquark fields which belong 
to (6,1,-2/3) representation of SU{3)c x SU{2)l x U{1)y- This along with the leptoquark 
in (3,2,1/6) representation makes the couplings blow up as they evolve at around 10^^ GeV, 
long before unification. This is true for all representation with higher dimensions. We have 
also investigated fiipped SU{^) x U{\) models pl where leptoquarks in 10 dimensional rep- 
resentation can arise from a product of 10 x e^. These states can give rise to leptoquarks 
in e~p scattering as well as e^p. Leptoquarks of the type Qe'^ and d'^e'^ are present. We 
will point out that their masses at the low energy scale are very close to each other in any 
realistic scenario. In this paper we give up grandunification for the most part, and instead 
consider the limited assumtion of gauge coupling unification at a large scale. 

We shall assume that there is a leptoquark field with transformation (3,2,1/6) and to 
remove anomalies, the conjugate representation is also present. The leptoquark field will 
couple to dc and (i/, e^), but the conjugate fields have no direct coupling to matter. The 
conjugate field has gauge couplings that will allow them to be produced in pp collisions at 
FERMILAB, but only indirectly through mixing at HERA. Since there are terms in the 
potential which allows mixing between the field and the conjugate field, the signature will 
bethe same as the original field. Consequently, if these masses are not well seperated, Teva- 
tron could have seen or will see the effect of this additional state in the near future. To 
maintain perturbative gauge unification we now add diquark and dilepton content of 10 
representation (3,1,-2/3) and (1,1,1) and their conjugate. We shall refer to these multiplets 
of new fields as / and their conjugate as /. Another interesting aspect of the R-parity 
conserving theory is the presence of the fermionic partners of the leptoquarks-the lepto- 
quarkino [^. The leptoquarkino decays into missing energy plus the same final state as a 
leptoquark does. The leptoquarkino production crosssection at Tevatron is higher than the 
leptoquarks for the same mass, since they are fermions. Consequently the leptoquarkino 
has to be heavier than the leptoquark to escape detection. This can be a problem since 
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our natuaral expectation is the fermion masses are smaller than the scalar masses. We will 
show that it is possible to have the fermion masses become larger than the scalar masses by 
adding a singlet field. Since we are not using the SU(5), we can couple the leptoquark to 
just the first generation, diquarks and dileptons need to couple to additional generations. 
We will consider the possibility that the leptoquark couples to other generation and derive 
bounds on the coupling to the matter fields. Our strategy is to write the superpotential 
involving the new fields and find the mass spectrum of both the fermions and scalars at low 
scale using appropriate renormalization group equations (RGEs) and the D-terms using two 
specific supersymmetric breaking models. We shall consider both the gauge mediated SUSY 
breaking scenario as well as gravity induced SUSY breaking since the process of transmission 
of information of SUSY breaking to the observable sector is still an area of active interest. 
We shall explore the experimental signature for the new fields in present and future colliders. 

We first show that with the choice of multiplets gauge unification is preserved even in 
two loop analysis. In Fig.l we show the two loop evolution of gauge couplngs assuming 
a — XjVn.^^ siT^Q-w = 0.2321 and = 0.118 at Mz- If either the diquark or dilepton 
is left out, the gauge unification is not possibile. As mentioned, a similar attempt with 
higher dimensional representation will cause coupling to diverge well below unification. In 
order to calculate the mass spectrum of this model we need a superpotential. Let us first 
discuss the superpotential in the gauge mediated symmetry breaking (GMSB) scenario. 
The gauge mediated SUSY breaking models have become very popular in last one year for 
several reasons. Since the SUSY breaking is communicated to the obsevable sector by SM 
gauge group or some other gauge group, the soft terms in these models are naturally fiavor 
symmeteric. These models also have less number of input parameters than the models where 
SUSY breaking is communicated by gravity eg. A (coefficent of the trilinear term) term and 
B term (coefficent of the bilinear term) in the potentials are at the GMSB sacle. The 
hghtest supersymmetric particle (LSP) in these models is gravitino and the next to LSP 
(NLSP) is either a stau or a neutralino, depending on the parameter space. The signals of 
these kinds of SUSY breaking are spectacular since they have hard photon (neutralino is 

4 



the NLSP) or a lot of taus (stau is the NLSP) in addition to the missing energy in the final 
states. 

The superpotential with the leptoquarks at the GMSB scale can be written as: 

M' 

+ —^1^1 

where 5*1 is a SM singlet, the necessity of the singlet field will be discussed later. Mj's are 
the soft breaking masses that get introduced by the same mechanism that generates the ix 
at the GMSB scale. We will make an ansatz that there is a common value for Mj when 
we perform the calculation. Similarly, to make the theory predictive, we assume the As 
are unified to a comman value of 1. When we write the potential from this superpotential 
there are terms like AjAj^j^jS'i and Mji?j0j0j. There are actually three more terms in the 
superpotential involving the leptoquark and the quark and lepton field, the diquark and the 
two quark fields and the dilepton and two lepton fields These terms arise from from a GUT 
scale term like hkLd'^ip. The coupling /i^required to explain the HERA events is quite small 
(~0.05). Since the diquarks are color triplet (antisymmetric combination of two color triplet 
fields) fields and the dilepton are SU(2) singlet field, they need at least two generations for 
this kind of coupling to survive. We can assume that these new field couple to first and the 
3rd generation. In fact if the leptoquarks {Ld'^) are coupled to the first generation only, the 
latest CDF bound rules out the leptoquark mass of about 210 GeV at 95% C.L. |T3] using 



the NLO QCD calculation |TJ]. If we assume that the leptoquark couples to all generations 
fiavor diagonally, then K ^ fie gets an unacceptable tree level contribution. We therefore 
forbid the leptoquark coupling to the second generation. However leptoquark can still couple 
to the third generation. This induces tree level contribution to -B ^ re. The bound we 
get is hi X < 5.7 x 10^'^, where hi and h^ are the leptoquark coupling to the first and 
the third generation. Consequently a third generation coupling to the leptoquark of the 
same magnitude as the first generation (~ 0.05) is allowed. Furthermore the inclusion of 
the third family implies that the leptoquark can decay into a r and a b quark. In this case 



the crossection needs to be multiplied by a factor of 1/4 and the leptoquark mass bound 
reduces to 190 GeV. 

There are also soft supersymmetry breaking gaugino (gluino, wino, bino) and scalar 
masses (squarks, sleptons, Higgs, fields in / and /) due the gauge mediated interaction with 
the messenger fields. These masses at the messenger scale M are given by [0,|n[ 



and 

m\M) = 2 H / (A) ^ a (^)' (3) 

where Oj, i = 1,2,3 are the three SM gauge couplings and ki = 1, 1,3/5 for SU(3), SU(2), 
and U(l), respectively. The Ci are zero for gauge singlets, and 4/3, 3/4, and (V/2)^ for 
the fundamental representations of SU{3) and SU{2) and U{1)y respectively. Here n cor- 
responds to number of (5 + 5) multiplets, and g{x) and f{x) are messenger scale threshold 
functions with x = A/M. The net scalar mass of the (pid'^ and (peui'Pd'^d'' fields have contribu- 
tion from Mi and as well as mg. The masses of the fields in / and / are same at the gauge 
mediated scale. We run these A's from the GUT scale to the GMSB scale by the following 
RGEs : 

2V\i = Ai(- ^ CliiATia.-:) + 5A^ + 3A^ + GA^), (4) 

i 

2VX2 = A2(- E C2i(47r«i) + + h\l + 6A^), 

i 

2VX3 = X3{-J2 C3i{47cai) + Xl + 3X1 + 8X1 



where i refers to the gauge group, V = ^^^4: and 



^ ' ' 5 / ' V 3 ' ' 15/ ' V 3 ' ' 15 



We need to run the soft masses, B terms, A terms from the GMSB scale down to the weak 
scale wfor the new fields as well. The RGB's for the soft masses are given as: 

/lf2 

Vm'^ = - E CU^^a,)M^ + X]{m] + m] + ^ + A]), (5) 



where j=l,2,3 represent Ldf^, d'^d'^, ve. We also have a similar equation for rrij. MiS are the 
guagino masses. The A's evolve according to the following RGE's: 

VALd^ = C3iiA7iai)Mi + GXIAlci^ + SA^Arfc^c + \lA,„ (6) 

i 

VAdcdc = Y.C2i{A7rai)Mi + SXjAid^ + SA^Arfc^c + XjA^^, 

i 

VA,e = E Cli{47iai)Mi + GXlALd^ + SXjA^.d^ + SXfA,, 

i 

Similarly, B's evolve according to the following RGE's: 

VBj = J2 Cj.(47ra,)M, + 2X]Aj (7) 

i 

The mass parameters M's evolve according to the RGBs: 

2VMj = M,{-J2 CMAnai) + 2A^') (8) 

i 

We can similarly write down the RGE for the soft breaking mass of the singlet 5*1 . We use 
the ref. [0] for the RGE's for the gauge couplings,for the other soft masses e.g the gaugino 
{Mi), squark and slepton masses and for the other parameters like At^^r- We keep fi as free 
parameter which is determined at the weak scale in the tree level by: 

Ml _Ml^-Ml^tan'^l3 



tari^fi — 1 



- (9) 



where tan/3 = ^. After we run these equations down to the weak scale we have to include 
the D term contribution to the scalar masses, which is of the form: 

Y F 3 

These D-terms cause a mass splitting between / and / scalars. Further, there is mixing 
between / and / scalars through the B term in the potential. This is a large effect for 
diquarks and leptoquarks but small for dileptons. This is because B term can be large for 
color fields at the weak scale even if its boundary value is at the GMSB scale. The actual 
mixing for the leptoquark is very nearly 45°. Consequently at HERA, the lighter combination 
of the mixed state is produced although e+g couples only to the f representation. The higher 



mass state is out of reach for HERA with the present energy. Similarly at Tevatron, one can 
produce either state through gluon exchange, but mass limits apply to to the lower mass 
state. The mass of the singlet field and the coupling of the singlet to the leptoquarks are 
needed here since if we look at RGE's we find the gauge part and the Yukawa part act in the 
opposite direction and we will want the leptoquark mass to be lower than the leptoquarkino 
mass. This is because the leptoquarkino production crosssection is always larger than the 
leptoquark production cross sections for the same mass, and the final state in the case of a 
leptoquarkino pair production will also have the leptons and jets. So the same bound which 
applies to leptoquark will also apply to the leptoquarkino. In Table 1 we show the masses 
of the new fields as well as the masses of the squarks, sleptons, neutralinos, charginos and 
gluino. The input parameters at the GMSB scale are A, M, n, tan/3, MjS and M' . The 
Yukawa couplings for the regular fermions are evaluated from their masses. The Yukawa 
couphng (Aj) for the leptoquarks are assumed to be unified at the GUT scale and its value 
at that scale is taken to be 1. At the GMSB scale Ai — 0.3, A2 = 0.7, A3 = 0.8. We have 
assumed the masses Mj (where i=Ld'^, d''d'^,a.iad ve) to be the same at the GMSB scale. We 
vary the common mass M and M' in such a way so as to get the leptoquark mass around 200- 
220 GeV. This fixes the diquark and dilepton masses also. Prom the table 1 it is easy to see 
that the leptoquark mass ~ 200 — 220 GeV can be obtained over a wide range of parameter 
space. In that parameter space the diquarks are much heavier than the leptoquarks and 
the leptoquarkinos are also heavier than the leptoquarks. The dileptons can be lighter or 
heavier. We can see from the table that the mass states of 0Lcfc are well seperated due to 
the large B term. This is also true for the diquarks but not for the dileptons where B is 
small. If we take the M^'s in the ratio of the Yukawa coupling (A^'s), the mass spectrum of 
the scalars would remain unchanged, only the dileptino mass will reduce to almost half the 
value showed in the tables. 

Let us discuss the signals of the new fields in different colliders. In Tevatron the diquarks 
will be produced pairwisc through s channel gluon exchange just like the squarks. Each 
diquark decays into a pair of jets. So we have 4 jets in the final state. The dileptons 
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can be seen at the Tevatron through Drell-Yan production with subsequent decay into 2 
leptons plus missing energy. Similar signals can be produced by the charged Higgs pair, 
W pair or the selectron pair. Leptoquark pair decays into 2jet+21epton. Beside Tevatron 
these particles can be pair produced in the electron-positron collider. One can also produce 
leptoquark and dilepton singly in e7 collider through a t channel exchange. The final 
state will be /+jet in the case of leptoquark and a 2 leptons plus missing energy (one lepton 
has to be other than electron, since the dileptons are connected to two leptons of different 
generations) in the case of dilepton. The fermionic partners of these fields: leptoquarkino, 
dileptino or diquarkino have interesting signatures in the colliders. The diquarkino will decay 
into a quark and a squark. The squark will eventually decay into a quark and neutralino 
(quark gluino channel is not suitable since gluino mass is large in GMSB models). The 
neutralino will give rise to a photon and a gravitino. When diquarkinos are pair produced 
the final state wil have 4 j+2 hard photons +missing energy. In the case of stau being the 
NLSP the diquarkino will decay into 4j+4r lepton +missing energy, since each neutralino 
will deacy 100% of the time into 2r +plus missing energy. The Prp distribution of those 4 
rs will be different and two of them will have very high Pt- This signature is very hard 
to miss. The production of these diquarkinos will be similar to that of top quark. The 
diquarkino mass of about 370 GeV will approximately give rise to one event (4jets + 4ror27 
plus missing energy) in Tevatron with a luminosity of 110 ph~^ . So far 6 events having 4 



jets plus missing energy (missing energy> 60 GeV have been observed [|T^, but they are 
claimed to be consistent with the SM precesses plus detector induced background. The 
dileptino will decay to a lepton and a sneutrino or neutrino and slepton and sneutrino or 
slepton will decay into neutrino or lepton and a neutralino, which will give rise to final states 
lepton+missingenergy (slepton mass is less than the neutralino mass) when stau is lightest 
or the lepton+7+ missing energy when neutralino is the lightest in the gauge mediated 
models. Consequently when they are produced in Tevatron or electron positron collider the 
dileptino pair can give rise to l^l^ + 27+missingenergy signal. One such event 9where both 
the leptons are electrons) was claimed to be observed recently fl^. The same signal can be 
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produced by the selections also. In an e7 collider the dileptino can be singly produced in 
the association of sneutrino through t-channel exchange of dileptino. Sneutrino decay gives 
rise to 7 plus missing energy (neutralino is NLSP) or 4r plus missing energy (lighter stau is 
the NLSP) in the final state. The leptoquarkino decays into a spositron (spositron decays 
into a positron and a neutralino with Xo ~^ iG) and d quark (or squark and an electron, 
but the squark masses are large in these models) with the final state in the leptoquarkino 
pair production process has e"'"e~77jj plus missing energy without any SM background. If 
lighter stau is the NLSP and neutralino is the NNLSP, the final state in the leptoquarkino 
pair production process has e+e^2r+2r^jj plus missing energy (since in this casexo 4r 
plus missing energy) . Out of these six leptons in the final state, one r pair (produced 
from the decay of stau) has much higher Pt than the other leptons. So far each CDF 
||T8| and DO |19| have reported one event in the e^e~]] plus missing energy channel. If we 



apply that bound we get leptoquarkino mass has to be greater than 345 GeV, since the pair 
producton crosssection is around 0.04 pb and e^e~ detection efficiency is 0.2 assuming the 
leptoquarkino branching ratio to be 1. If however leptoquarkino branching ratio is 1/2 the 
bound becomes 290 GeV. One could also look for the signals of leptoquarkino production 
in e+7 colliders. The leptoquarkino can be singly produced in such a collider. In the gauge 
mediated SUSY breaking scenario, the final state has either a hard photon or r+r^ along 
with electron +jets+missing energy. In all the above cases the fermion component of the 
new multiplet can be distinguished from the spin component by hard photon or excess rs 
and the missing energy. 

If we had used the leptoquarks in SU(5)xf/(l) unifying group, we would have 3 lep- 
toquarks and a dilepton as discussed before. In this case the above analysis would be 
unchanged. However all the leptoquarks will have almost the same soft symmetry breaking 
mass at the GMSB scale (since all of them have color quantum number) and their masses 
at the weak scale also will be very close (D term can not produce much of a splitting). This 
means that not only HERA, but also the Tevatron will see these 3 leptoquarks within a small 
mass range. The dilepton (z/^e^) can be produced at FERMILAB through Drell-Yan mech- 
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anism and has an interesting signature {e^ e" e'^ e" j j j j) arising from the decay of dilepton 
into ei^R, followed by vji eW~ — > ejj. 

In the supergravity motivated theory the soft masses are introduced at the GUT scale. 
We assume the squark, Higgs and slepton the new f multiplet all have same soft mass mo- 
The gaugino masses are also unified to a common mass mi/2- The other parameters are the 
bilinear coefficient of the // and the 5*1 5*1 term in the superpotential and the coefficient of 
the same bilnear terms in potential. Unlike the gauge mediated model the bilinear coefficent 
need not be to start with. In the case of supergravity motivated theory the signals are 
different from the GMSB models signals. Since neutralino is the LSP, there is no hard 
photon or high Pt ts in the final state. For example, the diquarkino would decay into quark 
and squark and the squark will deacy into a quark and missing energy. The final state will 
be 4jets+missing energy when the diquarkino is pair produced. The leptoquarkino pair will 
give rise to a jjl''"l~ plus missing energy and the dileptino pair will give rise to I'^'l" plus 
missing. The missing energy however would distinguish the fermion signal from the boson 
one. 

In Table 2 we show the masses of the new fields as well as the masses of the squarkss, 
sleptons, neutralinos, charginos and gluino. The input parameters at the GUT scale are 
the universal scalar mass mo, universal gaugino mass mi/2, Mj's, M', tan/?, A's and i?'s. 
Again for simplicity we have assumed the masses Mj to be same at the GUT scale. We 
also have assumed that other than the B associated with the f and /, all the S's are at 
the unification scale. We also assume that all the A's are at the GUT scale. From the 
table 2 it is easy to see that the leptoquark mass ~ 200 — 220 GeV can be obtained in wide 
range of parameter space. In that parameter space the diquarks are much heavier than the 
leptoquarks and the leptoquarkinos are also heavier than the leptoquarks. The dileptons 
can be lighter or heavier. The lighter combination of ^^^c and ^^^c is again much lighter 
than the heavier combination 

So far we have considered the case with just one type of leptoquark i.e If the 

charge current anamoly is also established, which is not at all clear from the present data, 
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we shall need an additional representation like 45 which has Lu'^ type of leptoquarks. These 
leptoquarks will then mix because of symmetry breaking with the Ld'^ type of leptoquarks 



and using the technique showed in the ref. one can get a charge current excess. If 
charged current events do hold up, such additional multiplets will need to be considered. 

In conclusion we have used leptoquark in the context of gauge unification and have 
discussed the mass spectrum and the signatures of the additionals fields needed to maintain 
the unification. We have showed that the fermionic partner of the leptoquark will be heavier 
than the scalar part. We have also showed that one combination of the leptoquark field and 
conjugate field (we need to add in order to cancel the anomaly) is lighter than the other 
combination, even though the soft mass generation mechanism is same for both types of 
fields. The heavier combination is still out of reach of Tevatron. 

We would like to thank K. S. Babu and S. Nandi for useful discussions. The work has 
been supported by a DOE grant no. DE-FG06-854ER-40224. 
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TABLE CAPTIONS 



Table 1 : Mass spectrum for the new fields (diquarks, leptoquarks and dileptons), stop and 
sbottom squarks {mi,m^), stau lepton(mf), chargino (x^)) lightest neutralino (m^o), 
gluino mass (m^) are shown in the gauge mediated supersymmetry breaking model. 
The scalar component of the new fields are written as mj and the fermionic components 
are written as M^. The two entries for the scalars give the hghter and the heavier mass 
eigenstates. 

Table 2 : Same mass spectrum as in Table 1 is shown for the supergravity scenario. 

FIGURE CAPTIONS 

Figure 1 : The evolution of the coupling constants are shown in presence of the new fields using 
two loop RGEs. . 
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Table 1 





Scenario 1 


Scenario2 


Scenario3 


Scenario4 


Scenario5 


masses 


A = 60 TeV, 


A = 60TeV, 


A = 60 TeV, 


A = 30 TeV, 


A = 25TeV, 


(GeV) 


n=l, M = 4A 


n=l, M = 4A 


n=l, M = lO^A 


n=2, M = lO^A 


n=2, M = lO^A 




tan /3=3 


tan /3=35 


tan j3=3 


tan/3=3 


tan/3=12 




371,580 


372,581 


346,630 


307,596 


271,579 




388 


388 


427 


418 


444 




152,174 


151,177 


207,225 


244,259 


270,287 


Me. 


305 


305 


309 


309 


311 




210,500 


210,510 


210,592 


210,570 


203,570 




210,500 


210,510 


210,592 


210,570 


203,570 




399 


400 


444 


434 


463 


m^o 


78 


80 


77 


75 


72 


m^± 


143,403 


148,347 


144,446 


135,346 


113,274 


nif 


111,222 


54,223 


119,230 


89,171 


65 


nit 


614,719 


633,705 


542,678 


441,564 


373,475 




685,726 


653,700 


638,655 


509,547 


414,435 




538 


538 


534 


535 


447 




-379 


-322 


-426 


-315 


-238 
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Table 2 





Scenario 1 


Scenario2 


Scenario3 


Scenario4 


Scenario5 


masses 

(GeV) 


mo = 250GeV, 
mi/2 = 250GeV, 
tan /3=3 


mo = 265GeV, 

mi/2=200GeV, 
tan /9=5 


mo = 280GeV, 
mi/2 = 180GeV, 
tan /3=30 


mo = 320GeV, 
mi/2 = 160GeV, 
tan/3=10 


mo = 230GeV, 
mi/2 = 240GeV 
tan/3=10 




350,789 


321,651 


304,605 


302,562 


333,727 




455 


401 


401 


402 


566 




183,306 


202,312 


216,324 


247,348 


119,261 


Me. 


204 


180 


180 


180 


202 




212,841 


212,692 


205,642 


204,593 


207,772 




212,841 


218,692 


205,642 


204,593 


207,772 




512 


401 


451 


452 


565 


m^o 


100 


80 


72 


64 


97 


m^± 


191,502 


148,389 


136,338 


322,345 


185,428 


nif 


269,311 


277,305 


248,316 


89,171 


244,293 


nit 


443,671 


372,567 


354,513 


329,501 


446,640 


mg 


596,650 


498,548 


429,498 


438,498 


570,619 




683 


547 


491 


437 


656 




-483 


-477 


-314 


-301 


-407 
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FIGURES 
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